Self-organization, scaling and collapse in a coupled 
automaton model of foragers and vegetation 
resources with seed dispersal 

D Boyer and O Lopez-Corona 

Institute de Fi'sica, Universidad Nacional Autonoma de Mexico, Apartado Postal 
20-364, 01000 Mexico D.F., Mexico 

E-mail: boyer@fisica.unam.mx 

Abstract. We introduce a model of traveling agents {e.g. frugivorous animals) who 
feed on randomly located vegetation patches and disperse their seeds, thus modifying 
the spatial distribution of resources in the long term. It is assumed that the survival 
probability of a seed increases with the distance to the parent patch and decreases with 
the size of the colonized patch. In turn, the foraging agents use a deterministic strategy 
with memory, that makes them visit the largest possible patches accessible within 
minimal travelling distances. The combination of these interactions produce complex 
spatio-temporal patterns. If the patches have a small initial size, the vegetation total 
mass (biomass) increases with time and reaches a maximum corresponding to a self- 
organized critical state with power-law distributed patch sizes and Levy-like movement 
patterns for the foragers. However, this state collapses as the biomass sharply decreases 
to reach a noisy stationary regime characterized by corrections to scaling. In systems 
with low plant competition, the efhciency of the foraging rules leads to the formation 
of heterogeneous vegetation patterns with 1//" frequency spectra, and contributes, 
rather counter-intuitively, to lower the biomass levels. 
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1. Introduction 

Animal movement and its ecological implications is a discipline that attracts a growing 
interest [1]. The study of animal displacements gives valuable clues on how complex 
organisms adapt to their environment, in particular to search, prepare and consume 
food Foraging problems have motivated many modeling approaches, sometimes 
inspired from the physics of the random walk problem [3]. In a way similar to 
anomalously diffusing particles in a physical context, the displacement patterns of a 
variety of animals (albatrosses bumble-bees [1], primates [5], gastropods [6], jackals 
[7], seals [8] and sharks [9], among others) involve many spatio-temporal scales and 
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are sometimes well described by Levy walks or intermittent processes with power-law 
distributions. On a biological point of view, wide fluctuations in the movements of 
a herbivorous or frugivorous animal are interesting as they may reflect a variety of 
behavioural responses induced by a complex environment with resources distributed 
heterogeneously [IDl [IH [El US] . 

Plant ecosystems are out-of-equilibrium and exhibit rich structures and dynamics. 
Spatial patterns in the distribution of plant species are highly non random and 
contain many characteristic length scales [H]. Patch distributions in rain forests 
have fractal properties, suggesting that these systems could be near a self-organized 
critical state driven by the slow growth of trees and sudden mortality avalanches |15j . 
Other observations report that tree sizes (and therefore fruit contents) in template 
and tropical forests are distributed according to inverse power-laws [T6l ITf] . When 
water resources are a limiting factor, continuous models show that plant interactions 
produce aggregation in patches that self-organize at larger scales to form more or less 
regular patterns [18], or disordered ones where patch areas obey scale-free probability 
distribution functions [T9] . 

Seed dispersal represents an important animal/plant interaction that may 
contribute to the formation of complex ecological patterns. Seed dispersal at long 
distances has been identified as an important structuring factor of tree communities 
[201 [2T] . Fruit eating animals [e.g., spider monkeys [22]) swallow the seeds of many tree 
species and deposit them, through faeces, practically intact and away from the parent 
tree after a transit time of a few hours. Between 60 and 90% of the seeds of tree species 
of tropical forests are dispersed by vertebrates that feed on fruit [23] , specially primates 
[211 [251 [26]. 

The aim of this article is to study an automaton model of moving foragers that 
modify, via seed dispersal, the long term structure of the resources they consume. In 
turn, these resources also determine the foragers displacements, who use cognitive skills 
to explore their medium in an efficient, non-random way. The model assumes that two 
factors influence the growth success of a seed: the distance to the parent plant [20] 
and competition due to the presence of other plants [21j. Despite that the model is 
over-simplifled, the forager /resources coupled dynamics leads to rich behaviours, like 
self-organized states with power-law statistics for patch sizes and animal movement 
lengths. 

The following two Sections describe the model and its background. The results and 
the underlying mechanism leading to power-laws in this system are presented in Section 
m and conclusions in Section [H 

2. Background: a foraging model without plant dynamics 

We flrst describe the movement rules of the model forager in a stationary distribution 
of resource patches (see refs. [131 l2Z])- Consider a two-dimensional square domain of 
area unity containing flxed, point-like patches randomly and uniformly distributed. 
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To each patch i is assigned a fixed size (or food content), ki = 1,2..., a integer drawn 
from a given distribution, p{k). 

Like many other animals, primates use cognitive maps to navigate their environment 
[28| [29] . Evidence shows that travels to fruiting trees are more frequent than as suggested 
by random null models Primates also keep record of the sites they have visited 

in a recent past [28]. For simplicity, we assume that our model forager has a perfect 
knowledge of the positions and sizes of all the patches in the system. Initially, a forager 
is located on a patch chosen at random. The following deterministic foraging rules are 
then iteratively applied at every time step (t ^ t + 1): 

(i) The forager located at patch i feeds on that patch, the fruit content 
decreasing by one unit: fcj — > fcj — 1. 

(a) If ki has reached the value 0, the forager chooses an other patch, j, such 
that kj/dij is maximal over all the allowed patches j i in the system, where kj is the 
food content of patch j and dij the Euchdean distance between patches i and j. With 
this rule, the next visited patch (the "best" patch) has a large food content and/or is 
at a short distance from i. We assume that the travel from i to j takes one time unit. 

(Hi) The forager does not revisit previously visited patches. 
This model produces complex trajectories that have been studied in details in refs. 
[T3] [27] and discussed in connection with spider monkeys foraging patterns observed in 
the field [5]. 

The model has a remarkable property, of interest in the following. Let us define 
the forager mean-displacement R(t) as (|r(t-|-to) ~i"(''^o)|) with r(t) the forager position 
at time t, the averages being taken over different times to and independent disorder 
realizations. At fixed patch number N and time t, if the resource size distribution is the 
inverse power-law p{k) = ck~^ with (3 = 3, then the mean displacement R{t) is maximal 
[T3] . In other words, media with this size distribution induce maximal displacements, 
see figure 1. (This property still holds if the forager travels at constant velocity instead 
of moving in one time unit from one patch to the other [13].) 

The feature above can be understood qualitatively by noting that if the medium 
is very homogeneous (say, (3^1), then all patches are similar in size: given rule 
(a), the forager chooses essentially nearby patches. Trajectories are thus composed of 
small steps and diffusion is relatively slow [13]. On the contrary, if the medium is very 
heterogeneous (/? ~ 1), patches with ki ^ 1 are numerous: the forager often performs a 
large step to reach a very good patch and stays there a long time feeding, given rule (z). 
The forager activity is dominated by these long trapping times, resulting in very slow 
(nearly frozen) diffusion. An intermediate situation corresponds to /3 = 3, for which 
the best patch from a given point is often far away (at distances much larger than the 
typical distance between nearest-neighbour patches), but these good patches still have 
reasonable sizes, so that the forager does not remain trapped feeding on them during 
very long periods of time. 

Let us note that, at the special resource exponent value /? = 3, the trajectory of the 
forager closely resembles a Levy flight. Numerical simulations show that the distribution 




Figure 1. Numerically obtained mean displacement (in unit oi Iq = N ^/^) as 
a function of time for the model forager in a random medium with resources size 
distributed as = Cfc-'^ {N = 10^). 



function of the distances separating successively visited patches is asymptotically given 
by the power-law [131 [27] 

p{i) ~ r'. (1) 

The foraging patterns of spider monkeys are well described by the distribution ([T]) [5lfT3]. 
Even more, these animals feed on trees whose size distribution obeys p{k) ~ with 
P ~ 2.6, a value close to 3 [13]. 



3. Modified model with plant dynamics 

In the previous model, forager motion is induced by the medium. In the generalization 
considered from now on, the forager follows the same rules but also modifies its 
environment through seed dispersal. Hence, the distribution of resource sizes, p{k), 
is no longer held fixed and can slowly change over time. We assume that foragers are 
the main mechanism of seed dispersal. 

At t = 0, N point-like patches, all with size fcj = 1, are randomly and uniformly 
distributed in the square domain of unit area. A forager initially located on a patch 
chosen at random follows the rules (i)-(iii) above. In addition, 

(iv) every time units (the digestion time), a seed is deposited at the patch 
where the forager is located; 

(v) every r^a«fc(^ Td) time units, the walk ends and the forager is removed; 
the patches are refreshed to their initial ki values; the patches that have received a seed 
that has survived (see below) increase their size by one unit, ki ^ ki + 1; 

(vi) a. new forager {i.e. not representing necessarily the same individual, but 
still having a perfect knowledge of the updated environment) is located on a patch 
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chosen at random and the process is iterated as above for another r^aik time units; 

(vii) a plant that has grown from a seed deposited time units earher 
(xm ^ Twalk) dics and the patch size decreases, ki ^ ki — 1; 

(via) in rule (vii), the size of a patch does not decrease below the minimal 
value ki = 1. 

We assume that two factors determine the survival probability of a deposited seed 
(or the growth success of a plant here) in stage (v). The first assumption is based 
on observations that seeds dispersed far away survive better, as they are less likely to 
attract seed predators and to be transmitted parasites or diseases from the parent plant 
[20] . The second assumption takes into account competition among plants for limited 
nutrients [21]. Let us note k the size of the patch where the seed is deposited and / its 
distance to the parent patch, which is the patch where the forager was located Td time 
units ago. The survival probability of the seed, Pg, is set to Pg = Pd{l)Pc{k) with 

p,{i) = [i/{i+io)r, (2) 

1 - iik<k 



10 if > kfnax) 

with Iq = N~^/'^ the typical distance between nearest-neighbour patches and k^ax a 
fixed integer accounting for competition effects. The parameter k^ax is proportional, 
say, to the nutrients concentration. According to ([2]), the survival probability increases 
from to 1 as the distance from the parent patch increases, whereas it decreases from 
1 to as the size of the colonized patch increases from 1 to kmax-, the maximum patch 
size. In the following, we fix n = 1 in Eq.(l2]), the results presented below being not 
qualitatively modified if larger values of n are chosen. 



4. Results 



It is natural to ask the following question. Given that in the model the survival 
probability increases with the distance to the parent plant and that the largest forager 
displacements are produced in media with the size distribution p{k) ~ /c~^, do resources 
self-organize towards this particular scaling-law (and foraging patterns towards the 
Levy law ([I]))? As a consequence of the memory-based foraging rules (i)-(iii), such 
environments should have the highest average biomass production rate. As the foragers 
start dispersing seeds, small heterogeneities are produced (fcj 7^ 1), and, from a 
dynamical system point of view, the fastest growing modes may dominate the dynamics 
asymptotically. 

Successful seeds contribute to the emergence of larger, more attractive patches that 
are also susceptible of being visited more often in the future. This positive feedback loop 
contributes to the increase of heterogeneities, within some limits. If very big patches 
are produced, due to rule (i) many seeds will be dispersed at the distance / = and will 
die, given the kernel ([2]). The other stabilizing effects are competition and mortality. 
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Figure 2. Time series (in units of r^aik) of the biomass, for two values of the 
competition parameter k^ax- The parameters are N = 10^, = 10, r^aik = 500 
and r,„ = 5 10"*. 

4.I. The rise and fall of scaling 

To investigate the possibility of the scenario sketched above, let us first consider the 
behaviour of the system biomass, M = YliLi ki, as a function of time. As displayed in 
figure [2], M initially increases and reaches a maximum value. It then suffers a abrupt 
drop, followed by a noisy stationary regime. The asymptotic average biomass depends 
strongly on plant competition: unexpectedly, it is lower at low competition levels (large 
)• 

The system actually builds spatial heterogeneities during the initial growth regime. 
Biomass is maximum at tmax = Tm = (lOOr^aifc, here), when the first plants grown from 
dispersed seeds start to die. As shown by figure [3], the size distribution p{k) at t = is 
perfectly fitted by the power-law k^^ for k <^ kmax, independently of the value of kmax- 
At the same time, and as expected from Eq.([l]), the forager step length distribution 
(figure H]) tends to the scaling-law P{1) ~ in a range of intermediate values of /, with 
a truncation at large / due to finite kmax, Tm and A^. Note the existence of steps of order 
of the system size (IOO/q, here). In some sense, the system self-organizes into a critical 
state of maximal dispersion. 

Interestingly, this scaling behaviour does not persist at larger times t ^ 
Tjn- Simultaneously to the biomass collapse (Figl2]), corrections to scaling appear. 
Asymptotically, the shape of the patch distribution decreases rapidly with patch size 
and exhibits a practically constant fat tail, see figure [3l Hence, most patches have a 
small characteristic size and coexist with a few "outliers" of size of order k^nr- This 
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Figure 3. Patch size distributions (in log-decimal representation) for kmax — 20 
(left) and kmax = 50 (right), at two different times: t^axi— lOOr^jaZfe, here), when the 
biomass is maximum (see figE]), and t = bOOOrwaik, in the stationary regime. Averages 
are performed over 10'^ independent initial conditions. The straight lines have a slope 
-3. 



1 
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Figure 4. Forager step length distributions at tmax, when biomass is maximal, for 
kmax — 20, 50 and 100. The straight line has a slope —2. 



separation into two characteristic sizes becomes more pronounced at large k^ax- 

Unexpectedly, dispersal by animals does not manage to stabilize the critical state 
for which it is responsible. Despite that the system rapidly self-organizes into a state 
with optimal seed survival, it becomes overpopulated: mortality is much higher than 
birth shortly after Tm and the size of many intermediate patches start to shrink. This 
feature is probably due to the fact that, at t = Tm, the system is much more crowded 
with plants than at t = 0. The birth rate is therefore lower than in the initial growth 
regime because of higher competition (Eq.(l3])): older dying plants are not replaced by 
the same quantity of new ones. In addition, due to their cognitive maps and the rules 
(i)-(iii), foragers neglect many small, unattractive patches, among them the shrinking 
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Figure 5. Left: Mean value of the biomass in the stationary regime for several values 
of the parameter kmax- Averages are performed over 100 independent runs. The curve 
has a maximum at kmax — 7 (corresponding to a fairly strong plant competition) and a 
bistable region centered around kmax = 58 (low plant competition). Same parameters 
as in Fig. El Right: Log-decimal power spectrum S{f) = (M(/)M(-/)), with M(/) 
the Fourier transform of the biomass M{t) in the stationary regime (see figure [2]), 
for kmax — 2, 7, 20 and 60. Averages are taken over 60 runs. The straight line is a 
power-law with exponent —1.4. 



ones, that would be suitable for plant growth. Instead, they keep revisiting a small 
number of nearly saturated large patches (outliers). 

4-2. Biomass in the stationary regime 

Asymptotically, foragers therefore concentrate their activity on a fraction of the available 
land, resulting in a reduced biomass. Figure [5] (left) displays the average biomass in the 
stationary regime as a function of the maximum allowed patch size, kmax (or nutrients 
concentration). Counter- intuitively, but in agreement with the comments above, the 
general tendency is a biomass decrease with increasing kmax- The biomass is large at 
small kmax (relatively low nutrients/high competition regime) and presents an extrema 
at kmax = 7. At this parameter value, patches are relatively homogeneous in size: via 
foragers visits, seeds can colonize many different patches. At higher kmax, foragers start 
to visit large patches preferentially, neglecting smaller ones. 

Supporting this interpretation, the power-spectra of the biomass time series M{t) 
reveal an increase in complexity with increasing kmax, see figure [5] (right). At high 
competition levels {kmax = 2), the fluctuations of M{t) around its mean value are due 
to births and deaths that occur in a roughly independent way, due to a homogeneous 
dispersion of seeds. The spectrum is that of a white noise. At low competition 
{kmax = 60), the power spectrum is well approximated by a power-law 1//" with 
a ~ 1.4 over more than one decade, indicating long-range temporal correlations. 
The spatially heterogeneous colonization of plants generates periods of high mortality 
( "avalanches" ) of widely varying durations, followed by periods of easier recolonization. 
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This dynamics is obviously reminiscent of the punctuated relaxation of sand-pile models 
in self-organized criticality |30j . 

An other unexpected phenomena is observed at low competition levels, in the 
interval 50 < kmax < 70 for the parameter values considered here. In this interval, the 
average biomass increases again and exhibits a second maximum at k^ax — 58. A closer 
look reveals that the system actually converges towards two distinct states, depending on 
the initial condition. At identical parameters values, the system sometimes ends up in a 
low, "base" biomass level that follows the tendency described above, and sometimes 
exhibits a significantly higher biomass ("excited" state). Biomass fluctuations are 
relatively small among different systems belonging to a same class (base or excited), 
which makes possible the computation of the mean values separately (see figure [5]). The 
origin of this bistability is unclear, although it is a known phenomenon in other models 
of vegetation pattern formation [18]. 

5. Discussion and conclusions 

We have proposed a new mechanism leading to the emergence of many spatio-temporal 
scales in the movement patterns of foraging animals. In the model proposed, foragers use 
mental maps to choose feeding patches and disperse seeds along their trajectories, thus 
affecting the long term distribution of their food resources. This simplified model focuses 
on plant-forager interactions, neglecting other important factors (wind, gravity...) of 
seed dispersal [21]. It is built on a generalization of a previous model, where scale- 
free displacement patterns of knowledgeable animals emerge from their interaction with 
resources that are distributed according to an a priori given power-law function p!3l[27] . 
In the field, the distributions of the movements of spider monkeys and of the size of their 
fruiting trees are in good agreement with those given by that model for a particular 
parameter value, where animal displacements are maxima [T3] . 

The present approach shows that a memory-based ranging behaviour generates 
highly heterogeneous seed deposition patterns, a conclusion also reached in ref. [26] 
with the use of a spatially explicit model parameterized with field-collected spider 
monkeys movement data. These findings suggest that, over large temporal scales, 
tree distributions may form complex spatial structures due to the presence of foraging 
animals. The present model proposes a theoretical test of this hypothesis: the 
distribution of resources is not held fixed and spatial heterogeneities self-organize 
spontaneously under the influence of positive feedback loops in the system dynamics. 

Other existing theories of Levy [1] or intermittent [311 [32] foraging assume that 
animals are memoryless and do not have any information on the location of resource 
patches. In these contexts, they execute a given Markovian stochastic processes to 
find preys (usually randomly distributed in space) that are detectable only at short 
distance. Movements with nontrivial distributions or rules are optimal for finding 
preys most efficiently in some cases. Whereas this approach can be justified for marine 
animals foraging in unpredictable environments [9|, frugivorous vertebrates rely on fixed 
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resources and memory plays an essential role Nevertheless, the introduction 

of limited knowledge and the use of search modes would improve the realism of the 
modelling approach presented here. 

Despite that, at each step, our model foragers maximize an efficiency function, their 
foraging activity can not be considered as optimal in the long term. The hypothesized 
relationship between long distance dispersal and species diversity in plant communities 
[13] is probably not a simple one [33J. The model exhibits sudden plant mortality 
avalanches that are consequences of a restricted land use and a lack of colonization 
of regions with low plant density. Counter-intuitively, biomass levels are lower when 
the conditions for colonization are favourable, i.e. at low plant competition levels. 
The same mechanism leading to the emergence of the self-organized critical state with 
optimal seed dispersal is also responsible for its rapid collapse. This ecological "crash" is 
the product of an intelligent foraging behaviour based on the satisfaction of immediate 
needs (feeding). Here, foragers do not change their strategy when plant resources start 
to shrink. It would be interesting to investigate whether the introduction of noise (or 
"irrationality") in the deterministic decision rules improves this situation. 

The model developed here is similar in some aspects to sand-pile models of self- 
organized criticality (SOC) [30]. The system is driven by successful dispersed seeds 
and biomass dissipated by mortality. SOC-like models have been proposed to explain 
large extinctions in the fossil record |34j and tree dynamics on ecological time-scales 
in rain forests [15]. In the latter example, l//-noise signals have been identified in 
simulated biomass time series [15] (unfortunately, biomass in real forests is very difficult 
to measure). Qualitatively similar results are obtained here with other assumptions. 
The stationary regime of our model, however, is not strictly speaking asymptotically 
critical since the shape of the patch size distribution contains two characteristic sizes: a 
small one, of order 1, and a much larger one, corresponding to the presence of "outliers". 
Such distributions seem to be ubiquitous in driven self-organized systems, though: they 
appear in earthquake models [35] and also in other models of plant dynamics in the 
presence of limited nutrients [19] . 

The robustness of the features described above should be tested by modifying the 
foraging rules and the dispersion kernels. A more realistic, resource-dependent forager 
demography is an aspect that should also be considered. 
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